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ABSTRACT 

We report on the discovery of a new heavily polluted white dwarf. The DAZ white 
dwarf GALEX J193156.8+011745 was identified in a joint GALEX/GSG survey of 
ultraviolet-excess objects. Optical spectra obtained at ESO NTT show strong ab- 
sorption lines of magnesium and silicon and a detailed abundance analysis based on 
VLT-Kueyen UVES spectra reveal super-solar abundances of silicon and magnesium, 
and near-solar abundances of oxygen, calcium, and iron. The overall abundance pat- 
tern bears the signature of ongoing accretion onto the white dwarf atmosphere. The 
infrared spectral energy distribution shows an excess in the H and K bands likely 
associated with the accretion source. 
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1 INTRODUCTION 

The star GALEX J193156. 8+01 1745 (hereafter 
GALEX J 193 1+0 11 7) is a hydrogen-rich white dwarf 
recently ident i fied in a joint ultraviolet/optical survey 
jVennes et all l2010h based on the Galaxy Evolution 
Explorer (GALEX) all-sky survejQ and the GSC2.3.2 cata- 
logue. The Third U.S. Naval Observatory CCD Astrograph 
Catalog locates the star at R.A.(2000)= 19 31 56.933, 
Dec.(2000)= +01 17 44.13 with a proper motion of 
fi a cos S — —77.6 ± 6.4 and [is = —6.1 ± 6.1 mas yr _1 . 
The distance modulus implies that the star is nearby 
(« 55 pc) and that it is imbedded in the Galactic plane 
(\z\ ~ 8 pc). Based on echelle spectroscopy we show that 
GALEX J1931+0117 is a DAZ white dwarf and is one of 
the most extreme cases of externally polluted white dwarf 
atmospheres. 

The DAZ white dwarf atmospheres are hydrogen- 
dominated with heavy element abundances ranging from 
nearly solar down to detection limits several orders of mag- 
nitud e below solar (|Zuckerman et all 120031 : iKoester et al.l 
l2005h . The DAZ and, more particularly, the DZ phenom- 
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ena are interpreted as evidence of on-going accretion of 
circumstellar material onto white dwarf atmospheres (e.g., 
iKilic fe Re dficld 200?]). The source of the material has been 



ascribed to comets (|Alcock. Fristrom. fc Siegelmanl 1 1986 
and more recently to tidally-disrupted asteroids (see lJural 
2008), although close, low-mass companions in post-c ommon 
envelope binaries are also a known source of material l|Debesl 
120061 : iKawka et al1l2008u . 

The observational evidence in favour of the accre- 
tion of debris material resides i n (1) the inferred com - 
position of the accretion flow (|Zuckerman et all l2007f ) 
showing an overabundance of refractor y elemen t s (se e 
a discussion on the solar system by iLodderj l2003l l. 
and (2) the direct spectroscopic signature of gaseous 
discs as in the cases of SDSS J122859.93+104032.9 
and SDSS J104341.53+085558.2 (rGansicke et all 120061 : 

or dusty discs (e .g. , 



Gansickc, Marsh, & Southworthl 1200 



Farihi. Jura, fc Z uckcrman 2009; Br inkworth et alj|200a ) 



The accreted material may also be supplied by mass 
loss from a close, low-mass, and possibly substellar com- 
panion. Only a few substellar companion to white dwarf 
stars are know n such as the whit e dwarf plus L8-9 resolved 
pair PHL 5028 (|Steele et al.ll2009h or clo se DA plus L8 post- 
common envelope bi nary WD0137-349 l|Maxted et al.ll2006l : 
iBurleigh et al.l2006lh Relatively luminous DAZ wh ite dwarfs 
such as EG 102 (|Holberg. Barstow. fc Greenlll997n may hide 
low-m ass companions, but iDebes. Sigurdsso n. fc Woodgatel 
(2005) limits them to substellar types. The mass loss ma- 
terial captured by the deep potential well of a white dwarf 
should reflect the composition of the chromosphere of the 
red dwarf, which, in most cases, should be close to solar. 
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Table 1. Photometry 
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Figure 1. Model fits (solid lines) to the low dispersion (EFOSC2) 
Balmer line profiles (jagged lines) and best-fit parameter values 
(T^, log fl). 



In this context, the discovery of the peculiar DAZ white 
dwarf GALEX J1931+0117 is timely. We present a series of 
ultraviolet (UV), optical and near infrared (NIR) observa- 
tions (Section 2) and our model atmosphere analysis (Sec- 
tion 3.1) revealing a DAZ white dwarf with a peculiar abun- 
dance pattern (Section 3.2). A notable NIR excess (Section 
3.3) may hold clues concerning the nature of this external 
source. We summarise and conclude in Section 4. 



2 OBSERVATIONS 

IVennes et al.l (|2010| ) describe a selection of hot sub-luminous 
stars in a joint GALEX/GSC2.3.2 survey. Briefly, a UV ex- 
cess along with measurable proper-motion is interpreted as 
a likely, but not definitive, characteristic of hot subdwarfs 
(sdB, sdO) or white dwarfs. To confirm the selection, follow- 
up spectroscopy was done for a number of candidates includ- 
ing GALEX J1931+0117. We obtained two low-dispersion 
spectra (t cxp = 600 s) using EFOSC2 attached to the New 
Technology Telescope (NTT) at La Silla Observatory on UT 
2009 August 25 (grating #11, dispersion = 4.15 A pixel" 1 ). 
The spectral coverage is from 3680 to 7408A at a resolution 
of w 15.8A. 

Following the identification of GALEX J1931+0117 as a 
DAZ white dwarf we obtained three echelle spectra (i cxp = 
1450 s) using UVES attached to the VLT-Kueyen on UT 
2009 November 12 (dispersion ~0.02 A pixel" 1 ). The first 
two echelle spectra were obtained with the red arm and cover 
the ranges 5698-7532A and 7660-9464A. A third spectrum 
obtained with the blue arm covers the range 3757-4985A. 

We collected the JHK p hotometry from the 2MASS 
survey (|Skrutskie et alj 120061 ) and the GALEX photom - 
etry, corrected for nonlinearity (|Morrissev et al.l l2007h . 
from the GR4 ±5 release. Despite a reddening in dex of 
E B -v = 0.4054 (|Schlegel. Finkbeiner. fe Davidll998h in the 



Bandpass 


Measurement 


GALEX FUV 


13.20±0.10 


GALEX NUV 


13.55±0.10 


V 


14.20±0.02 


2MASS J 


14.66±0.05 


2MASS H 


14.55±0.09 


2MASS K 


14.45±0.10 



line of sight toward GALEX J1931+0117, the star is rel- 
atively nearby and does not show the effect of reddening 



(Ee 



< 0.01). Finally, we used the acquisition images 



(V filter) at NTT to measure the apparent magnitude of 
GALEX J1931+0117. We employed the acquisition images 
of the flux calibration standard Feige 110 (V = 11.832) and 
EG 21 (V = 11.40) obtained on the same night at the NTT 
to calibrate the measured counts and derive the V magni- 
tudes. Table [1] summarises the available photometric mea- 
surements. The field surrounding GALEX J1931+0117 is 
somewhat crowded with objects 4.8 and 8.3" away. These 
separations are sufficient to exclude significant flux contam- 
ination in the optical and infrared images. The ultraviolet 
(FUV and NUV) images are possibly contaminated. How- 
ever, the optical/infrared indices imply that the nearby stars 
are not significant ultraviolet sources. 



3 ANALYSIS AND DISCUSSION 
3.1 Atmospheric parameters 

First, we constrain the effective temperature (T e ff) and the 
surface gravity (logg). Two separate measurements show 
a slight systematic difference in log g. Fitting the low- 
dispersion Ha to H9 profiles we measure T e s — 20880 ± 240 
and logg = 7.84 ±0.04, while fitting the high-dispersion H/3 
to H8 we measure T cff = 20890 ±140 and log 5 = 7.93 ±0.03. 
Table [2] lists the weighted averages. The error bars for the 
surface gravity were extended in order to encompass both of 
our measurements. We computed the corresponding mass, 
age, an d absolute V magnitude usin g the mass-radius rela- 
tions of lBenvenuto fc Althausl <| 19991 ) with a carbon/oxygen 
core, thick hydrogen layer QogMn/M* = —4), and zero 
metallicity. Finally we estimated the distance using the pho- 
tometric parallax method. Figure Q] shows the analysis of the 
low-dispersion Balmer line profiles. 

Next, we calculated the space motion of the star using 
the proper-motion (Third U.S. Naval Observatory CCD As- 
trograph Catalog) and the barycentric velocity (see Table[3]) 
corrected for the gravitational redshift (j B — 25.9^24)) 
«bary,corr — 12.5 km s -1 . We empl o yed t he methodology de- 



scribed b y I Johnson fc Soderbloml ( 19871 ) and the solar mo- 
tion from iHogg et al l (|2005h . The star belongs to the thin 
disc population. 

Our basic model atmosphere grid was computed in lo- 
cal thermodynamic equilibrium (LTE), but additional non- 
LTE model atmospheres were computed usin g TLUSTY 
l|Hubenv fc Land 1 19951 ; lLanz fc Hubert Il995h . We deter- 
mined the non-LTE metallicity vector by fitting a non-LTE 
model at 20900 K and log g — 7.9 to a grid of pure-hydrogen 
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Table 2. Properties of the DAZ GALEX J1931+0117 



Parameter 



Measurement 



Effective temperature 
Surface gravity 
Mass 

Absolute magnitude 
Cooling age 

Distance 
Kinematics 



: 20890 ± 120 K 

n +0.03 



T c ft 

log ff =7.90lg;g| (c.g.s.) 
M = 0.57±°'^ M 
M V = 10.50+J!'^ 
i cool = 4.5 ± 0.6 X 10 7 yrs 



d = 55 
(f, V, IV) 



+ 3 



(28,5,22) km s" 



Table 3. UVES line identification 



Ion 


A 


E.W. 


^bary 




(A) 


(mA) 


(km s — 1 ) 


Sill 


3853.660 


17. 


36.6 


Sill 


3856.020 


56. 


36.8 


Sill 


3862.609 


50. 


36.7 


Can 


3933.663 


50. 


37.1 


Can 


3968.469 


11. 


37.7 


Sin 


4128.070 


56. 


36.4 


Sin 


4130.890 


83. 


36.2 


Mgii 


4481. 13/. 33 


579. 


33.7 


Fen 


4549.474 


6. 


35.3 


Sim 


4552.622 


11. 


38.1 


Fen 


4923.927 


12. 


37.6 


Sin 


5957.560 


34. 


40.8 


Sin 


5978.930 


19. 


43.2 


Sin 


6347.100 


132. 


36.4 


Sin 


6371.360 


102. 


36.4 


Oi 


7771.944 


39. 


39.2 


Oi 


7774.166 


31. 


39.0 


Oi 


7775.388 


22. 


41.5 


Mgn 


7877.054 


46. 


45.5 


Mgii 


7896.366 


103. 


45.5 


Can 


8542.091 


18. 


35.8 



a Average barycentric velocity t>bary = 38.4 ±3.2 km a -1 . 

LTE models. We measured AT off = +120 K and A log 3 = 
—0.04. We conclude that despite the high-metallicity of the 
atmosphere, the systematic error introduced by using pure- 
H LTE models is within statistical errors. 



3.2 Photospheric abundances and diffusion 

Table |3] lists the heavy-element lines identified in the UVES 
spectra. The average barycentric velocity of the trace ele- 
ments («bary = 38.4 ± 3.2 km s -1 ) is in agreement with the 
Ha velocity = 37.6 ± 0.8 kms -1 . We determined the 
photospheric abundances by fitting the observed line pro- 
files with non-LTE models. Figure [2] shows selected heavy 
element line profiles and best-fit models. We did not notice 
line splitting due to a magnetic field at the resolution limit 
(R = A/AA « 40,000). Since AA/A = 4.7 x 10 -7 AB where 
A is in A and B in MG, the limit implied is B < 10 kG. 

Table [4] lists the measured abundances. Helium and 
carbon are not detected and we employed HelA5875.7 and 
C IIA4267. 001-183 to set abundance upper limits. The in- 
frared calcium triplet region shows weak photospheric ab- 
sorption but no gaseous disc emission with a limit of E.W. 
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Figure 2. Echclle spectra (UVES) of the DAZ 
GALEX J1931+0117 showing a strong Mg doublet and 
weaker oxygen, silicon, calcium and iron lines. 

lA. The abundance of magnesium in GALEX J1931+0117 
exceeds that of SDSS J122859.93+104032.9 by a factor of 
5 but no evidence of a gaseous disc is found in disagree- 
ment with a correlation between magnesium abundance and 
cal cium emission strength in hot DAZ wh ite dwarfs f ound 
by IGansicke. Marsh, fc Southworthl l|2007| ). However, ljural 
(2008) argues that gaseous discs are rare and that in t he two 
cases discussed by IGansicke. Marsh. &: Southworthl (I2007T ) 
an additional source of heating must be responsible for this 
unusual phenomenon. 

The near-solar heavy element abundance in the atmo- 
sphere of GALEX J1931+0117 must be supplied from an 
external source. The accreted material reaches optical depth 
unity, tr » 1, and, assuming near steady-state, appears in 
the optical spectra at the measured abundances, n(X)/n(H). 
The elements are supplied in the accretion flow in pro- 
portions n(X)/n(H)fl ow that are not necessarily solar. The 
steady-state abundances are proportional to (1) the supplied 
abundances and (2) their respective diffusion time-scales. 
For example, the measured abundance of element X relative 
to silicon is given by 



n(X) _ r(X) n(X) 



(1) 



n(Si) r(Si) n(Si) flow 

where the time-scale r(X) is given at a selected point in the 
atmosphere by 
m 
p v 



r(X) 



(2) 



where m is the column mass density (gem -2 ), p is the lo- 
cal mass density (gem -3 ), and v is the diffusion velocity 
(cms - ) at that point. These quantities are calculated at 
tr = 1. Both m and p are provided by the model atmo- 
sphere. The diffusion veloci ty v = v(T,n p ,Xj) is calcul ated 
following the prescription of lFontaine fc Michaudl (1 19791 ) in a 
fully ionized hydrogen gas (Zn = 1), and the remaining vari- 
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ables, temperature T, proton density n p , and the ionization 
fractions Xi of element X are also provided by the model 
atmosphere. The diffusi on velocity is averaged over a ll ion- 
ization species following iMontmerle fc Michaudl §976). The 
relatively high concentrations of neutral helium and oxygen 
contribute toward shorter diffusion time scales compared to 
elements with lower concentrations of neutral species. 

Table [4] summarises our calculations. Using Equation 
(1) we calculated the abundances relative to silicon that 
are required in the accretion flow to reproduce the observed 
abundance ratios and compared them to solar ratios. 

Again, assuming steady state accretion, the total mass 
of element X accreted per year and mixed into an atmo- 
sphere of mass M 
radius, is given by 



AttR m, where R is the white dwarf 



M = A(X) 



n(X) AtvR 2 



(3) 



'n(H) r(X) 

where A(X) is the atomic weight of element X and the 
abundance is given relative to the main constituent. For 
example, magnesium and silicon must be accreted at a 
rate of 4 x fO 15 and 6 x fO 15 gyr" 1 , respectively. The to- 
tal accretion rate for magnesium and silicon is well within 
the range p redicted by the tidal destruction of asteroids 
ij Jural [200a ). but the relatively low abundance of calcium 
argues against this scenario. Abundance meas urements in 
the d usty helium- rich white dw arfs GD 40 (|Klein et al.l 
l20ld ) and GD 3 62 jKoesteill2009l ). and hydrogen-rich white 
dwarf G 29-38 (|Koesteri 120091 ) imply a calcium-rich accre- 
tion flow with n(Ca)/n(Mg) ~ O.f — 1. The accretion flow 
in GALEX J1931+0117 is low in calcium (n(Ca)/n(Mg) ~ 
0.02), somewhat below the solar ratio. Within uncertainties 
the accretion flow has near solar composition. 

Alternatively, assuming that the accreted material is 
supplied by mass loss from a low-mass companion, the ac- 
creted mass of the main constituent (hydrogen) is scaled on 
the accreted mass of trace elements. Averaging the observed 
variations between elements, the total accretion rate is esti- 
mated at Mace = 3 x 10 -14 Mq yr _1 . Using the post-common 
envelope white dwarf plus brown dwarf binary WD0137— 349 
as an example we arbitraril y locate a L5 c ompanion (see 
section 3.3) of 0.05 M ( see lKirkpatrickll2005t ) at a separa- 
tion of a = 1 Rq . The fraction of material captured by the 
white dwarf would then dictate the required mass-loss rate 
for the companion. Assuming in turn a geometrical cross- 
section, a gravitationally-focused cross-section, or fluid ac - 
cretion jWesemael. Henry, fc Shipmanl 1 19841 ; iDebed l2006h . 
the fraction is estimated at 5 x 10" 5 , 3 x 10" 3 , or 0.02-0.08, 
respectively. Under the assumption of fluid accretion, the 
required mass-loss rate is M acc = 0.4— 1.5 x 10 -12 Mq yr _1 . 
This rate exceeds the solar rate or the upper limit derived 
for Proxima Centauri implying that mass loss may be stim- 
ulated by tid al or irradiat ion effects in a close binary (see a 
discussion bv lDebesll2006r i. 



3.3 Infrared excess: low-mass star or dusty disc 

Figure [3] compares the spectral energy distribution from the 
UV to the NIR to computed composite spectra. Adopting 
class averages of Mj = 1 2.0, 13.0, 13.5, an d 14.0 for L2, L4, 
L5, and L6, respectively (Kirkpatrick 2005T). we scaled corre- 
sponding L templates from the IAC M-L-T Dwarf Catalog 
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Figure 3. UV to NIR spectral energy distribution based on 
GALEX FUV and NUV photometry, the NTT V magnitude, and 
2MASS JHK photometry compared to a white dwarf model (full 
line). Contributions from (top to bottom dashed lines) a L2, L4, 
L5, and L6 dwarf are added to the white dwarf model. 



ijMartfn. Cabrera, fc Cenizd [2005) to a distance of 55 pc, 
and added the flux to the projected NIR continuum of the 
DAZ white dwarf. The measured excess would corr espond 
to a L5 dwarf with T cff ~ 1700 K l|Kirkpatrickll2005tl . 

Assuming M2 = 0.05 Mq and a = IRq, the predicted 
L5 and DAZ velocity semi-amplitudes are = 314 ■ sini 
kms -1 and Koaz = 28 ■ sini kms -1 , respectively. Radial 
velocity shifts in the white dwarf spectra may be detectable 
in future observations. The presence of a close L dwarf com- 
panion should lead to notable optical/infrared light varia- 
tions. In particular, variable Ha emission is a noted spec- 
troscopic signa ture of post-common envelope systems (e.g., 
WD0137-349, IMaxted et "ail l2006t ) . but it is not apparent 
in any of our spectra. 

We cannot exclude the presence of a dusty circumstellar 
environment in addition to a L dwarf companion. However, 
by showing a steeper decline over the H- and K-bands the 
NIR spectral energy distribution of GALEX J1931+0117 
appears different than that of the warm, dusty whit e 
dwarf SDSS J122859.93+104032.9 (|Brinkworth et al.ll2009l ). 
If present, dust should reside at a distance from the star 
where the equilibrium temperature between dust and stellar 
radiation is below the sublimation temperature, T su b (see 
Jurall2008T l. For T sub between 1350 and 1650 K (see lLoddersl 
20031) the minimum distance from GALEX J1931+0117 or 
SDSS J 122859.93+104032 9 is b etween 1.7 and 1.1 R Q . 
In fact, iBrinkworth et all |2009) found that the gaseous 
and dusty discs in SDSS J122859. 93+104032.9 co-exist 
at a distance < 1.2 Rq from the white dwarf. Although 
GALEX J1931+0117 and SDSS J122859.93+104032.9 have 
nearly identical luminosities, the lack of a gaseous disc in 
the former and the dissimilarity of their NIR spectral en- 
ergy distributions imply different environments. 
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Table 4. Photosphcric abundances 



X 


n(X)/n(H) 


n(X)/n(Si) 


n(X)/n(Si) a 


r(Si)/r(X) » 


n(X)/n(Si) How c 


Departure from solar d 


He 


< 5 x icr 4 


< 7.8 


2.6 x 10 3 


60. 


< 5 x 10 2 


< 0.2 


C 


< 7 x icr 5 


< 1.1 


7.6 


1.07 


S 1-2 


$ 0.2 





2.4 ±0.2 X 1CT 4 


3.8 ±0.4 


14. 


7.4 


~ 30 


~ 2 


Mg 


7.1 ±0.3 x 10" 5 


1.1 ± 0.1 


1.1 


0.81 


- 0.9 


~ 0.8 


Si 


6.4 ±0.3 x 10~ 5 












Ca 


9.0 ±1.8 x 10~ 7 


0.014 ±0.003 


0.063 


1.23 


~ 0.02 


~ 0.3 


Fe 


2.9 ±0.6 x 10~ 5 


0.45 ± 0.09 


0.87 


1.86 


- 0.8 


~ 1 



a Abundances in the solar photosphere from Grcvcssc . Asplund. & Sau val (2007). b Ratio of diffusion time scales at tr = 1. 
c Predicted abundance ratio in the accretion flow. d Predicted abundance ratio in the accretion flow relative to solar. 



4 SUMMARY AND CONCLUSIONS 

We identified and analysed the properties of one of the most 
heavily polluted white dwarfs known. The source of the ma- 
terial remains unknown, although the measured NIR excess 
may be attributed to a L5 dwarf or to a warm debris disc. 
The measured abundances suggest that the material is ac- 
creted onto the atmosphere in solar proportions and favour 
a model involving a L dwarf companion in a close orbit. 
However, the Ha emission notable in such systems is absent 
in GALEX J1931±0117 and the mass loss imposed on the 
companion appears excessive. NIR intermediate dispersion 
spectroscopy should help determine the nature of the ac- 
cretion source. Additional optical spectra will be useful to 
trace putative orbital motions while Space Telescope Imag- 
ing Spectrograph high-dispersion spectra will be useful to 
constrain further the carbon abundance and extend the pat- 
tern to less abundant elements. 
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